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The coupling of spin and valley physics is nowadays regarded as a promising route toward next- 
generation spintronic and valleytronic devices. In the aim of engineering functional properties for 
valleytronic applications, we focus on the ferroelectric heterostructure BiAlOs/BilrOs, where the 
complex interplay among trigonal crystal field, layer degrees of freedom and spin-orbit coupling 
mediates a strong spin-valley coupling. Furthermore, we show that ferroelectricity provides a non¬ 
volatile handle to manipulate and switch the emerging valley-contrasting spin polarization. 
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Technological advances in future electronics strongly 
rely on the investigation of quantum degrees of freedom 
of electrons. Beside electron charge, spin represents the 
most studied example, due to its obvious connection with 
magnetic information storage. Spin may be also actively 
manipulated in devices based on spin-polarized transport 
[l[, as in the prototypical spin-FET proposed by Datta 
and DasQ- Relativistic spin-orbit coupling (SOC), pro¬ 
viding a link between the electron spin and orbital an¬ 
gular momenta, plays a central role in spin-polarization 
effects while preserving time-reversal symmetry. Such 
spin-polarization effects may appear even in centrosym- 
metric crystals [1, |^, provided that the ions experiencing 
a strong SOC have an inversion asymmetric local environ¬ 
ment, ultimately due to the localized character of SOC it¬ 
self. In the absence of bulk inversion symmetry, SOC lifts 
the spin degeneracy, leading to spin-ralitting phenomena 
known as DresselhausQ and Rashbap effects. In the lat¬ 
ter case, typically realized at interfaces or surfaces where 
inversion symmetry is structurally broken, the spin po¬ 
larization is induced by an effective fc-dependent mag¬ 
netic field, whose strength can be modulated by apply¬ 
ing an external electric field @,0- Alternatively, it has 
been recently understood that bulk Rashba-like effects 
can be realized in ferroelectric (FE) semiconductors 
additionally providing a non-volatile functionality (fer¬ 
roelectricity) which is tightly bound to spin polarization 
and, as such, allows for its full-electric control^. 

Other degrees of freedom have been recently addressed, 
most notably valley pseudospin, labelling the degenerate 
energy extrema in the momentum space, which could be 
used in valleytronic devices exploiting the valley index 
of carriers to process information [l(ij| . The research in 
this field has been boosted by the discovery of hexago¬ 
nal 2D crystals [lll[l2j|. such as graphene and transition- 
metal dichalcogenide monolayers, displaying valleys at 
the corners (K points) of the Brillouin zone. Various 
schemes have been proposed to generate vallw currents 
in graphene, using line defect s 11311 or strain [1^, and in 
biased dichalcogenide bilayers [ISj. The detailed investi¬ 


gation of these materials has significantly contributed to 
the fundamental understanding of valley physics, leading 
to the identification of intrins ic phys ical properties asso¬ 
ciated with valley occupancy [la-Ujj. Interestingly, the 
symmetry properties of valley pseudospin implies that 
valley-contrasting physical quantities may emerge when¬ 
ever inversion symmetry is broken [ibL Il7j | . In this re¬ 
spect, the possibility to permanently control the valley 
phenomenology via a non-volatile FE polarization would 
open the path to novel device paradigms, in analogy with 
the proposed FE Rashba semiconductors^. 

Most of the well-known FE materials are transition- 
metal oxides with perovskite structure. As recently 
pointed out, the hexagonal geometry can indeed be real¬ 
ized in perovskite heterostructures (AB 03 )„/(AB’ 03)2 
grown along the [111] direction, where the B’ ions in 
the (AB’ 03)2 bilayer sit on a corrugated graphene-like 
honeycomb lattice [H, [l^, as shown in Fig. [TJ Re¬ 
cent technological developments [2l| allow now to prepare 
such artificial layered heterostructures with atomic-scale 
precision [ 2 ^, offering an unprecedented playground for 
the realization of oxide-based advanced electronics. In 
fact, these (111) superlattices could allow the integration 
of functional properties typically found in oxides (23| with 
the spin-valley physics, which is expected to emerge in 
the hexagonal layered structure. Specifically, the “layer 
degree of freedom”, naturally emerging in the bilayer ge¬ 
ometry, can be described as a “pseudospin” up (down) 
labeling the state where the charge carrier is located in 
the upper (lower) layer, being therefore related to electri¬ 
cal polarization. Such a layer pseudospin is expected to 
couple with spin and valley pseudospin, possibly giving 
rise to exotic magnetoelectric effects [2^. 

Since the spin-valley contrasting properties can emerge 
when inversion symmetry is broken, in this letter we fo¬ 
cus on the quest for FE spin-valley physics in perovskite 
oxides. As design criteria, a robust FE system with polar¬ 
ization perpendicular to the hexagonal bilayer is needed. 
Since the BiFeOs-type rhombohedral crystal structure, 
displaying a [111] polar axis, is considered as an ideal 
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candidate to fulfill this criterion, we choose the nonmag¬ 
netic isostructural BiAlOa 23 recently synthesized, with 
a measured polarization Pg « 27 /rC/cm^ and a high crit¬ 
ical temperature, > 520° [ 2 ^, Aiming at a sizable 
spin splitting in a nonmagnetic system, in this study we 
propose the low-spin Ir^+ (5d®) as a candidate B’ ion, 
due to its large atomic SOC. Density-Functional-Theory 
(DFT) calculations were performed using the VASP code 
with GGA-PBE potential [ 2 ^. 



FIG. 1: (i) Multilayer structure of (BiIr 03 ) 2 (BiA 103 ) 4 . Oc¬ 
tahedral tilting and polarization direction are shown by ar¬ 
rows. (ii) Ir-O bond lengths in the bilayered structure and 
(iii) top view of the corrugated honeycomb Ir lattice. Elec¬ 
tronic bandstructure of FE (BiIr 03 ) 2 (BiA 103)4 without SOC 
(a) and with SOC (b). With SOC, P is tuned to P = 0 in a 
PE parent structure (c) and to P = —Pa (d), where Ps is the 
calculated spontaneous polarization. The Sz polarization is 
highlighted by blue (up) and red (down) color. Spin-splitting 
at VBM is indicated by block arrows. 

We first briefly summarize the outcome of our calcula¬ 
tions for hypothetical bulk BilrOs. We found it to sta¬ 
bilize in the polar R3c structure, with a calculated po¬ 
larization Pg = 70.6 fiC/cm^. As compared to the ideal 
perovskite structure, BilrOs displays two main distor¬ 
tion modes, i.e., IrOe octahedral tilting and Bi-0 polar 
distortion, the latter being responsible for the onset of 
ferroelectricity through the Bi lone-pair mechanism. 
alongside a compressive trigonal distortion. Ir^+-5(i® 
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FIG. 2: (a) Spin valleys of (BiIr 03 ) 2 (BiA 103 ) 4 : ±Sz polariza¬ 
tion projected on E{k) curve of the VBM band. Top of Pvbm 
is set as 0 eV. (b) FE polarization Pz and spin polarization Sz 
of the VBM band at K point as a function of the polar struc¬ 
tural distortion A. (c), (d) Spin texture of the VBM band for 
-|-Ps and -Ps crystal structure, respectively: Sx. and Sy com¬ 
ponents are shown by arrows, while Sz component is plotted 
by blue (-l-Sz) and red (—Sz) color of VBM band state. The 
first Brillouin zone is shown as a hexagon. 

electrons stabilize in a low-spin state, whereas the trig¬ 
onal distortion removes the degeneracy of fully-occupied 
t 2 g levels, which are split in a non-degenerate aig and 
twofold degenerate states. 

The bilayer structure displays a polar P3 symmetry, 
where Ir ions sit on a corrugated honeycomb lattice con¬ 
sisting of two trigonal layers intercalated by a single O 
layer, as shown in Figsllji) and (iii). With respect to 
bulk BilrOa, an additional structural distortion devel¬ 
ops, consisting in different distances of Ir layers from the 
central and the interfacial O layers, the latter showing 
on average longer Ir-O bonds [see FigjTKii)]. The par¬ 
ent paraelectric (PE) structure, belonging to the P321 
acentric (albeit non-polar, due to D 3 = 2 C 3 3 C 2 sym¬ 

metry) space group, is reached by tuning the offcenter- 
ing of the central O layer with respect to the Ir ones, 
other distortional modes showing minor changes. There- 
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fore, spin-splitting effects may appear also in the acentric 
PE bilayer. Due to the C 3 site-symmetry of Ir ions in 
both FE and PE structures, a built-in bulk Rashba-like 
(denoted as R-1) and Dresselhaus-like (D-l) effect can 
be anticipated, as arising from the local dipole field and 
acentric environment, respectively Q; in the presence of 
polar offcentering, additional FE-modulated R-l and D- 
1 spin polarizations should appear on top of such built-in 
effects. 

The electronic band-structure in the valence manifold 
mostly arises from the upper and lower Ir-(aig, e'g) or¬ 
bital states, forming molecular orbitals (MOs). Figures[T] 
(a)-(d) show the two bunches of d electron bands, namely 
the bonding and anti-bonding MOs arising from the d or¬ 
bitals at two Ir sites. All the bands are separated due to 
both trigonal crystal field (CF) and SOC-induced spin¬ 
splitting. The valence band maximum (VBM) is at the 
K point; remarkably, all valence bands at K are spin- 
polarized, displaying a large = 0.41 polarization 
at VBM, whereas and Sy components are zero. Due 
to time-reversal symmetry, relating K and K' points, 
the VBM shows a positive (negative) spin-polarization 
at the K {K') point, a manifestation of spin-valley cou¬ 
pling emerging due to the polar distortion, as shown 
in Figl2l(a). Interestingly, such valley-dependent spin- 
polarization appears to strongly depend on the FE phase, 
being switched when the FE polarization of the het¬ 
erostructure is reversed [Figs. [D (b) and (d)]. The size 
of the valley-spin polarization evolves proportionally to 
the size of the Bi-0 offcenter, i.e. of P, as shown in 
FigI21[b), vanishing when P = 0. In fact, even though a 
strong spin-splitting is still visible in the PE band struc¬ 
ture shown in FiglT] (c), a finite Sz spin polarization is 
found only along the M-T line, being forbidden by the 
C2 symmetries along the T-K line and at the valleys K, 
K' . However, a hidden valley-spin polarization can be 
identified when projecting the total Sz onto Ir ions be¬ 
longing to different layers, hinting to an intrinsic strong 
spin-valley coupling. 

Eventually, a complex pattern of in-plane spin polar¬ 
ization is shown in Figl^ (c) and (d), arising from the 
coexistence of several acentric distortions. The in-plane 
spin textures display the same chirality around the val¬ 
leys, suggesting a substantial decoupling from the valley 
pseudospin. When the FE polarization is switched, a 
complete reversal of the spin-polarization texture is not 
realized, even though a clear switching of the chirality 
can be seen in Figl2] (c) and (d); this suggests an inter¬ 
play of FE-dependent and built-in i?-l and D-l effects, 
the latter depending on the distortional modes which are 
almost unchanged when the polar offcentering is tuned 
between opposite Ps and persisting in the PE struc¬ 
ture. In summary, our ab-initio calculations show a com¬ 
plex coexistence of valley-dependent and bulk Rashba- 
like (valley-independent) spin-polarization effects in bi- 
layered BilrOa; remarkably, the spin-valley physics ap¬ 


pears to be tightly bound to the FE distortion, emerging 
only in the polar structures. 

The microscopic origin of the valley-spin polarization 
and its interplay with the FE polarization can be qualita¬ 
tively grasped by considering a simplified tight-binding 
model H = Pl^op + Hd -f Ffsoc0, describing t 2 g elec¬ 
trons moving on a corrugated honeycomb lattice in the 
presence of both SOC and trigonal CF. We will con¬ 
sider here only the effects of the polar distortion — in 
terms of the central O layer offcentering Sz — that af¬ 
fects both Ir-0 bond lengths and directions and thus 
determines the nearest-neighbor oxygen-mediated hop¬ 
ping interactions, neglecting the more complex struc¬ 
tural distortions unveiled by our realistic DFT calcula¬ 
tions. The Hamiltonian can be conveniently recast in a 
symmetry-adapted basis where the trigonal CF term be¬ 
comes diagonal, while the on-site SOC displays the well- 
known expression found for p (^ = 1) orbitals, the aig {e'g) 
states transforming as Pz {Px, Py) onesji^. Within this 
formulation, the analogies with graphene-like p-electron 
materials clearly emerge, the band-splitting induced by 
the trigonal CF mimicking the sp^ orbitals with cr, tt 
character, whereas SOC may lead to Rashba-like inter¬ 
actions in the corrugated honeycomb lattice due to spin¬ 
mixing terms coupling the tt (oig) state with the a (e'g) 
manifold. This is in striking contrast with M 0 S 2 
and d-electron dichalcogenide monolayers, where spin-up 
and spin-down states remain completely decoupled even 
in the presence of SOC due to the different trigonal pris¬ 
matic CF acting on d orbitals (see Table |I|. On the other 
hand, oxygens in the oxide bilayer provide an additional 
way by which inversion symmetry in the graphene-like 
honeycomb lattice can be broken. 

TABLE I: Bulk symmetry (space group SG), atomic-site sym¬ 
metry (site point group SPG), atomic SOG and expected spin- 
polarization effects in coupled spin-valley 2D compounds. Fol¬ 
lowing Ref. [^, R-l and D-l denote Rashba and Dresselhaus 
effects, while R-2 and D-2 denote the hidden spin polarizations 
due to the bulk IS. The electronic basis used for the SOC term 
is c = {px,Py,Pz), comprising cr, tt bands from sp^ valence 
bonds, for p-electron systems, while it is c = (eJ^j, e(g), ai(g)) 
for d-electron M 0 S 2 (BilrOs bilayer) as determined by the 
trigonal prismatic (octahedral-|-trigonal) CF. 
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A standard downfolding procedure onto the aig states 
near the K{K') point allows to derive an effective Hamil¬ 
tonian for the valence bands, iL®® = pSzPll^^ where 

the polar distortion has been included up to linear order 
in Sz [ 2 ^. Here, describes the low-energy electronic 
properties of the ideal (non-polar) perovskite bilayer, be- 
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ing 


i/o® = t±{Tko,Sa: + kySy) + Ta^SJo 

~k ^nij^xky o'ykx^Sz^ ( 1 ) 

where r = ±1 is a valley index labeling K and K', and cr 
and S denote Pauli matrices for the real spin and layer 
pseudospin, respectively, while t± , and Xr are effective 
parameters depending on hopping amplitudes, atomic 
SOC and trigonal Cf[ 23. As anticipated before, Eq. 
0 coincides with the low-energy model that has been 
derived for 2D atomic crystals with trigonal symmetry 
as silicene[3^. The second term describes the well-known 
spin-valley-layer coupling which is responsible for 

the opening of a gap in the valence manifold. The last 
term describes a layer-dependent Rashba-like coupling 
which arises from the SOC-induced mixing of spin-up 
and spin-down states mediated by inter-orbital hoppings 
in the corrugated honeycomb lattice. Due to the cen¬ 
tric symmetry of the ideal perovskite bilayer, no net spin 
polarization may appear; however, both the spin-valley- 
layer and the layer-Rashba couplings give rise to hidden 
layer-dependent spin polarizations In the presence of the 
polar distortion modeled by the central O layer offcen- 
tering, the low-energy electronic properties are further 
described by 

Sz + B^TCT;, a^TSxCTy + SyCTx) 

“t" {{SxCTx rSyCF y')ky -f {SyCJx -f T Sx^T y ) ] 

T XRi^Cxky eTykx^-i (2) 


where the additional effective parameters Ej,, a, A/ji 
can be expressed in terms of the tight-binding ones The 
coupling between spin and layer pseudospin described by 
i/®®, removing all the degeneracies in the valence man¬ 
ifold, can be loosely regarded at as a valley-dependent 
“magnetoelectric” interaction 2^ , explaining the emer¬ 
gence of a net spin polarization. It is also clear that, as 
highlighted by the linear dependence on the polar offcen- 
tering 5z, all coupling terms in Eq.([2]) change sign upon 
reversal of the EE polarization, causing a global switching 
of the spin polarization. Specifically, the averaged total 
spin-polarization can be derived at valleys K (K’), giving 
Sz = {o-z) = tSz{Ez + {E^ bI)^ (4a)2. De¬ 

spite the minimal parametrization of the starting tight- 
binding model, the derived effective Hamiltonian clearly 
shows that the EE distortion is responsible for both a 
valley-spin polarization and an i?-l effect, both modulated 
by the size of the polar distortion. 

In conclusion, we explored the possibility of engineer¬ 
ing the spin-valley physics in EE transition-metal ox¬ 
ide heterostructures. We have explicitly shown that 
the interplay of trigonal CE effects and SOC in (111) 
oxide bilayers leads to graphene-like low-energy elec¬ 
tronic properties, thus allowing for the appearance of 
coupled spin-valley physics. Interestingly, the spin- 


valley-layer coupling term can be regarded as mediat¬ 
ing an effective valley-dependent “magnetoelectric” in¬ 
teraction which manifests in complex patterns of spin 
and layer polarizations; when the polar distortion is in¬ 
duced by the oxygen-layer offcentering, a large valley 
spin-polarization, tuned by the EE polarization, devel¬ 
ops, as we verified for a proposed candidate heterostruc¬ 
ture comprising a robust EE BiAlOa and a bilayered 
BilrOa. The appearance of valley-contrasting properties 
and the “magnetoelectric” coupling brought in by the 
complex interplay of valley, spin and layer pseudospin 
have been proposed to have important implications in 
valleytronic devices|ld.ll7l.l24135|. Even though the pro¬ 
posed material might be not for the optical and electrical 
generation of spin and valley polarized carriers, due to its 
indirect band gap and the nearly flat dispersion of the 
valence band maximum, we believe that our theoretical 
findings are more general, suggesting that the realization 
of spin-valley physics in oxide heterostructures is indeed 
possible. Our theoretical analysis further suggests that 
such realization would allow in principle for larger effects 
(due to atomic SOC, which is typically large in 4d or 5d 
transition-metal ions), increased tunability (brought in 
by oxygens, both determining the trigonal CF splittings 
and mediating the hopping interactions) and for the inte¬ 
gration of additional functionalities, such as ferroelectric- 
ity, which could be exploited in advanced next-generation 
electronic devices. 
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